Submitted to ApJ 



SMALL STRUCTURES VIA THERMAL INSTABILITY OF 
PARTIALLY IONIZED PLASMA I. CONDENSATION MODE 

TSUBASA FUKUE AND HIDEYUKI KAMAYA ^ 
Department of Astronomy, School of Science, Kyoto University, Kyoto 606-8502, Japan 

tsubasaOkusastro .kyoto-u. ac . jp 
ABSTRACT 

Thermal instability of partially ionized plasma is investigated by means of a 
linear perturbation analysis. According to the previous studies under the one 
fluid approach, the thermal instability is suppressed due to the magnetic pres- 
sure. However, the previous studies did not precisely consider the effect of the 
ion-neutral friction, since they did not treat the flow as two fluid which is com- 
posed of ions and neutrals. Then, we revisit the effect of the ion-neutral friction 
of the two fluid to the growth of the thermal instability. According to our study, 
the characteristic features of the instability are the following four points: (1) 
The instability which is characterized by the mean molecular weight of neutrals 
is suppressed via the ion-neutral friction only when the magnetic field and the 
friction are sufficiently strong. The suppression owing to the friction occurs even 
along the field line. If the magnetic field and the friction are not so strong, the 
instability is not stabilized. (2) The effect of the friction and the magnetic field 
is mainly reduction of the growth rate of the thermal instability of weakly ion- 
ized plasma. (3) The effect of friction does not affect the critical wavelength Ap 
for the thermal instability. This yields that Ap of the weakly ionized plasma is 
not enlarged even when the magnetic field exists. We insist that the thermal 
instability of the weakly ionized plasma in the magnetic field can grow up even 
at the small length scale where the instability under the assumption of the one 
fluid plasma can not grow owing to the stabilization by the magnetic field. (4) 
The wavelength of the maximum growth rate of the instability shifts shortward 
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according to the decrement of the growth rate. This is because the friction is 
effective at rather larger scale. Therefore, smaller structures are expected to ap- 
pear than those without the ion-neutral friction. Our results indicate the friction 
with the magnetic field affects the morphology and evolution of the interstellar 
matter. In summary, the ion-neutral friction is important for the evaluation of 
the thermal instability in weakly ionized plasma along and perpendicular to the 
magnetic field. 

Subject headings: instabilities — ISM: evolution 



1. INTRODUCTION 



The recent progress of the interstellar medium (ISM) observations has established that 
the small and tiny scale structures of ISM is very ubiquitous. Historically, the initial obser- 
vational target is the neutral phase of ISM. Indeed, the tiny scale structures a re discovered 



by 21 cm absorption line observations against quasars with VLBI techniques (jPieter et al. 



19761 ). The results are confirmed by [Diamond et al.l (119891 ). To look for the small scale 
structures in the cold neutral phase of the ISM furtherm ore, the observat io ns of 21 cm 



absorp tion lines have been also performed against pulsars (IFrail et al.l 1 19941 ). iMeyer et al. 



( 119961 ) observe optical spectral lines to find small structures against clos e binary stars. Th e 
images of the small-scale H I have been taken with the MERLIN array (jPavis et al.lll996l ). 



These results are established by higher angular resolution observations with VLBA and VLA 
toward quasars (IFaison et al.lll998l : iFaison fc Gossll200ll ). which show small clumps on the 
order of a few AU i n neutral Galactic H I clouds. Co ld H I clouds have significant structure 
in subparsec scales ( ICibsonl 120001 : iBrogan et al.ll2005l ). Such small and tiny scale structures 
has been detected in the loc al interstellar medium a s H I absorption lines, although their 
column density is very small (iBraun &: Kanekarl 120051 ). 



In addition to the small and tiny scale structures of H I clo ud, the picture of the 



planetary nebula NGC 7293 shows fine small structures and knots ( iRodriguez et al.l 12002 



O'Dell et al.ll2004r). Even in the s tarforming regions, there are variety small scale structures. 
For example, iLanger et al.l (119951 ) observe some clumps with size from 0.007 to 0.021 pc in 
the Taurus Molecular Cloud 1 (TMCl), in particular, core D. The mass of these fragments is 
estimated to be < 0.01 — 0.15 M(t^. SPITZ ER has begun producing higher spatial resolution 
mid-infrared maps (jChurchwell et al. and revealed the fine structure of the starforming 

region. We think that it is possible for proto-brown dwarfs of < 0.08 Mq to exist there. A 
high resolution observation in future is expected to reveal hidden, small structures, as well 
as substellar objects with very small mass. In any ways, the tiny-scale structure seems to be 
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ubiquitous, not associated with large extinction (IHeilesI 119971 ) . 



To study the origin of the small-scale structures of ISM is important to understand 
the evolution of and structure formation in ISM. Especially, in the starforming regions, the 
small, low-mass structures can relate to the low-mass cutoff of the initial mass function, 
the coagulation unit to form the massive stars, and so on. Then, we should investigate the 
physical origin of these small and tiny structures in the partially ionized medium, since the 
cold H I and molecular clouds are weakly ionized. By the way, it is not easy for the small and 
tiny scale structure to form as a result of gra vitational ins t ability since their size is much 
smaller than the Jeans length. According to iLanger et al.l (Il995l ). indeed, the small-scale 
structures appear to be gravitationally unbound. This suggests that some fragmentation 
mechanisms but pure Jeans gravitational instability may be important in the clouds. We 
expect the thermal instability as this mechanism, and revisit it in this paper. 



The basics of the thermal instability has been summarized in iFieldl ( 1l965l ). When the 
following condition is satisfied, the system is thermally unstable: If the cooling becomes 
efficient as the temperature decreases, the thermal energy gets lost more and more. If the 
cooling becomes efficient during the fluid contraction, the system shrinks on and on. Impor- 
tantly, the critical length scale of the thermal instability is smaller than that of the dynamical 
instability like the Jeans instability. That is, even if a system is stable against the gravita- 
tional instability, the system can become thermally unstable. Then, the thermal instability 
can b e the physical origin o f the smaller-scale formation than the dynamical instability. In- 
deed, iKoyama fc Inutsuka propose that the dumpiness in clouds emerges naturally 
from their formation through the thermal instability. Their two-dimensional calculations 
follow the fragmentation into small cloudlets that result from the thermal instability in a 
shock-compressed layer. iBurkert fc Lirj (120001 ) investigate the cooling and fragmentation of 
optically thin gas with a power-law cooling function. According to them, small-scale pertur- 
bations have the potential to reach higher amplitude than large-scale fluctuations. Thermal 
instability can be important for the structure formation. 

We sketch ISM structure formation such as molecular clouds as follows. The hot ionized 
ISM cools to be cold and weakly ionized clouds. Ionization degree of ISM decreases as 
it evolves, and then the ISM fluid is often partially ionized. In partially ionized fluid, 
an ion component and neutral component interact each other, exchanging the momentum 
(i.e., ion- neutral friction). Especially, the partially ionized plasma in a magnetic fleld could 
not be treated as one fluid. This is because although the ion component in the partially 
ionized plasma is directly influenced by a magnetic fleld, the neutral component does not 
directly feel a magnetic fleld. In weakly ionized fluid, the neutral component is affected 
by a magnetic fleld via the ion-neutral friction. For example, ambipolar diffusion takes an 
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important role in dynamical evolution ow' i ng to the gravitational instability of ISM into 
protostar (IMestel &: Spitzerl Il956l : iNakand Il976l ). If the fluid frozen in a magnetic field 
contracts without ambipolar diffusion, the magnetic field bec omes too large and suppresses 
the growth by the magnetic pressure and tension. In addition, iNakand (119791 ) points out the 
possibility that ISM with sufficient magnetic flux quasistatically evolves into protostar. He 
also notices that the time scale of plasma drift depends on the ionization degree. The relation 
among the ion-neutral friction, the magnetic field, and the amount of the ion component is 
important to understand the detailed process. 

We here emphasize that the MHD approximation is not always applicable for our purpose 
in this paper. The thermal instability is effective in small structure formation, during which 
a neutral component is not always enough frozen in an ion component dynamically to be 
treated as one fluid. There is a possibility that the ion component can drift away owing to 
ambipolar diffusion, and that the system can contract owing to thermal instability. Then, 
in the present paper, to study how the ion-neutral drag and magnetic field influence the 
growth of the thermal instability, we treat the plasma as the two fluid of ions and neutrals. 
We assume that the ionization degree is very small, because we are interested in the final 
stage of the formation of a molecular cloud, in which the ionization degree decreases very 
much. We focus on the condensation mode of thermal instability, since it is important in 
structure formation, rather than the oscillation (overstable) mode. In §2, the problem is 
formulated. The property of the dispersion relation of the instability is presented in §3. 
Several discussions relating to the instability are found in §4. Applicability of our results to 
cold H I medium is briefly examined in §5, and then we summarize the paper in §6. 



2. BASIC EQUATIONS AND DISPERSION RELATION 



2.1. Basic Equations 

The basic equations for ion-neutral two component fluid are as follows. The continuity 
equation, equation of motion, energy equation, and equation of state are 



dpn 
dt 



+ V(pnf n) 



7-1 



dt 



Pn 

7 Pn 
7 - IPn 



dt 

'dpn 

dt 



+ (Vn ■ V)pn 



0, (1) 

-Vpn - Pnl^niiVn - Vi), (2) 
-p„An + V(i^nVTn), (3) 
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and 



n 

/in' 



for a neutral component. Those for an ion component are 
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dt 



+ V(pi'i;i) 



(4) 
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and the induction equation is 



dB 

'dt 



V X (vi X B). 



(9) 



Here, the subscript, n, denotes neutral, and subscript, i, does ion. Variables p, v, p, T, and 
-B express density, velocity, pressure, temperature, and a magnetic field, respectively. The 
cooling function A is defined as energy losses minus energy gains per gram per second. In 
addition, the coefficient of thermal conductivity is expressed by the gas constant by 7^, 
and the mean molecular weight by /i. 

The second terms on the right-hand sides of eq uations (fZ]) and ([6]) m ean the effect of 



dragging through the neutral- ion collision, similarly to lWatson et al.l (120041 ). The neutral- ion 
collision frequency i^ni is expressed as 7i/Pi, where 7,^ is the colhsion rate coefficient per unit 
mass, 

7. - (10) 

m-i -|- mn 

Here, a is the collision cross-section, v is the relative velocity between ion and neutral com- 
ponents, m is mass, and the angular bracket of () expresses a mean value over all velocities. 
It is pointed out that the larger the friction is, the better the one fluid approximation of the 
partially ionized plasma is like a classical MHD approximation. 

The effect of the heat flow between the two components due to the difference of the 
temperatures is neglected as second order terms because we assume the temperatures of two 
components under the equilibrium are identical. In addition, the frictional heating is also 
neglected as second order terms because the mutual velocity between the two components 
is zero in the equilibrium. 
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2.2. Dispersion Relation 

We solve the perturbed basic equations, considering the first order. The unperturbed 
state is supposed to be infinite, uniform, static, and isothermal in equilibrium. It is described 
with subscript 0, such as = p^o, Pn = Pno, = v^o = 0, = T^o, pi = Pio, Pi = Pio, 
Vi = = 0, Tj = Tio, and B = Bq. The magnetic field Bq in equilibrium can be taken 
as Bq = (0, 0, Bq) for simplicity. The net transfer of thermal energy for each component in 
equilibrium is zero, such as An(pno,^no) = and Ai(pio,Tio) = 0. 

Someone might think our unperturbed state is inadequate since the size of the unper- 
turbed medium is infinite. However, in our so-called WKB framework, we can state how 
our unperturbed condition is useful. As long as the WKB approximation is adopted, the 
whole of the system should be large. This can correspond to the fact that the scale length of 
the galactic gas disk is about 100 pc, which is always much longer than the so-called Field 
length, Ap. That is, since the interesting wave-length of the instability for the structure 
formation is smaller than the thickness of the galactic disk, our study is meaningful. 

We take the perturbation as 

a{r,t) = aiexp{nt + ik ■ r), (11) 



according to iFieldl (119651 ). where ai is the amplitude of perturbation, n is the growth rate 
of perturbation, and k is the wave number of perturbation. When n is real and positive, 
the perturbation grows exponentially as a condensation mode. If n is real and negative, 
the perturbation damps. When n is complex number, the system is oscillatory growing or 
damping. Of course, when ^[n] is positive, the system is so-called overstable. 

In order to know how the magnetic field affects on weakly ionized fluid, we mainly 
investigate the condensation mode, because we are interested in the formation and evolution 
of molecular clouds. We assume that 



Bq ± k and k \\ \\ v^i (12) 
for comfort. Due to the induction equation these assumptions lead to 

Bq II Bi, (13) 

[see equation (1A9P in Appendix |A] . It is emphasized that the assumed magnetic fleld is 
perpendicular to the fluid motion. 

Here, it is noticed that we can examine some other modes, which do not satisfy the 
condition (IT^ and (fT5]) . when Bq = 0. This means Bi = 0, owing to the induction equation 
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([9]). In other words, we can study the property of the dispersion relation when the magnetic 
field does not work effectively. This is interesting point in our analysis since those modes 
with Bq = yield understanding of the pure effect by the ion-neutral friction. It is also 
noted that there is no thermal instability in a pure Alfven mode which satisfies f ± fc, 
because the Alfven mode is effectively incompressible, and there is no net thermal exchange, 
even if the fluid is compressible. 

In this paper, since we are interested in the structure formation from H I medium, the 
fluid is assumed to be weakly ionized, i.e., pio ^ Pno- We also assume that temperatures of 
neutral and ion components are the same in the equilibrium, i.e., Tno = Tjo, because they are 
well mixed and collide enough to have the same temperature for fewness of ion component. 
Using these conditions and replacing \k\ = k, here is the dispersion relation (see Appendix 
|X]for the derivation in detail): 
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where 
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and they are evaluated in the equihbrium state. 

The dispersion relation (HM is the equation of sixth order of n. This dispersion relation 
connotes the two sets of three modes. One of the three modes in each set is condensation 
mode, and the other two modes are oscillation modes originally, if there is no friction effect. 
The first bracket on the left-hand side in equation (fT^ is the dispersion relation of three 
modes of only neutral fluid (IBip . The second in equation (IMI) is the dispersion relation of only 
ionized fluid in the magnetic field (lB2p . The right-hand side of equation fll4p expresses the 
momentum exchange through the ion-neutral drag. The properties of the thermal instability 
of one fluid approximation is also briefly mentioned in Appendix [Bl 



3. RESULTS 

3.1. Critical Wavelength of the Condensation Mode 

Observing the dispersion relation (fT4|) . we obtain the first result that a critical wave- 
length for the instability in a condensation mode is not affected by the friction. The critical 
wavelength is evaluated if we let n be zero, because the growth rate of n is zero at the critical 
wavelength. Thus, letting = 0, we find the right-hand side of equation f|T^ . which denotes 
the effect of the ion-neutral friction, diminishes. That is, the effect of the friction does not 
play any roles to determine the critical wavelength of the condensation mode. The crite- 
rion of the instability for the partially ionized plasma is estimated as if there is no friction. 
Although the critical wavelength is not affected, the friction affects the growth rate of the 
thermal instability. It is noted that the criterion of the overstability modes is not derived by 
the same way like the condensation mode since n is imaginary. 



3.2. Growth Rate of Condensation Mode of the Instability 

3.2.1. Description of the Conditions 

The dispersion relation f[T^ is numerically solved by means of MATHEMATICA, after 
the normalization (see Appendix [C]) . We study the condensation mode which is important 
on structure formation. To present the numerical dispersion relations, we take the following 
parameters : 7 = 5/3, an = ^nx/^np = 1/ 2, a-, = k\T/ k\p = 1/2, /9n = k^p/k^K = 0.01, 



Pi = kip/kiK = 0.01, for the comparison to iFieldl (119651 1. We set these parameters to be 



fixed because we are concerned for the basic nature of two-component fluid and compare 
the results with one-component fluid. We assume that /in = 1, /Xi = 1/2, for e.g., only H 
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atom plasma, x = 0.01, and k^p = k^p. Under these conditions, the dispersion relation ([T] 
is evaluated for the various quantities of v^/vsn which corresponds to the strength of the 
magnetic field, and z^nio/(^np^^sn) which corresponds to the strength of the friction. 

Before showing the precise results of the numerical dispersion relations, we examine the 
physical meaning of the axis of the following figures. In each figure, the horizontal axis is 
the wave number k in units of k^p. The vertical axis corresponds to the growth rate 3?[n] 
in units of k^pV^n- These normalizations are adopted to denote the characteristic property 
of the thermal instability. The wavelength is taken as A. We shall start considering the 
physical meaning of k^p. According to equation f|T6|) with approximation of PnoAnp ~ An, 

= /^n(7 - iKoA., 

7^^;sn^no ^ ^ 

PnoAn 1 
• 

PnO "^sn 

where Pno/PnoAn expresses the cooling time. Then, we find k~p is interpreted as the distance 
for the fluid to acoustically move within the cooling time. In other words, it is interpreted as 
the length scale in which the fluid can respond to the perturbation. Therefore, regarding the 
horizontal axis, {k/k^p)~^ ~ Afcnp is interpreted as the scalelength in units where the fluid 
can respond to the thermal instability within the cooling time (i.e., the growing time). We 
get some physical intuition once the above normalization is adopted. For example, we get 
physical intuition when {k/k^p)~^ ^ 1, the fluid is hard to respond to the thermal instability. 
This is because the system is too large and the acoustic wave which carries the disturbance 
can not travel through the whole system within the cooling time. Then, the growth rate 
at (k/k^p)^^ ^ 1 is reduced down in each figure. When {k/k^p)~^ < 1, the fluid is easy to 
respond to the thermal instability. 

Similarly, regarding the vertical axis, ^[n]/k^pVsn is interpreted as n multiplied by the 
cooling time. In other word, it is interpreted as the cooling time divided by the time scale 
for the growth of the instability. Then, the vertical axis means how fast the instability can 
grow up compared with the cooling time. The thermal instability is due to the cooling, then 
the value of ^[n]/k^pVsn should be less than unity as presented in each figure. 



3.2.2. Various Conditions of the Friction and Magnetic Field 

Figure [1] displays the dispersion relation of the real part of n when there exists neither 
the magnetic field nor friction. A dot-dashed curve corresponds to a mode characterized 
by Pi and a dashed one to that of Pn. No friction means that the motion of the ion and 
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neutral components are independent completely, and thus a dot-dashed curve corresponds 
to a single ionized fluid and a dashed one to a simple neutral fluid. The d ifferen c e bet ween 



the two curves comes from the mean molecular weight. This result follows iFieldl (Il965l ) 



Figure [2] shows the growth rate when there is no magnetic field but some friction [i.e., 
va/vsh = 0, i^nio/ {k^pVsn) = 0.03]. This condition yields understanding of the pure effect of 
the friction. When there exists the friction, a pure mode of an ion component or a neutral 
one does not exist. Therefore, interpretation is required on what a curve in the figure means 
when there exists the friction. The resultant modes with the friction are displays as two 
thick lines. For comparison, the two thin lines when there exists neither a magnetic field 
nor friction in Figure [1] are also plotted. A thick dot-dashed curve originates from an ion 
component mode of Figure 1. A thick dashed curve originates from an neutral component of 
Figure 1. These interpretations are reasonable since the effect of mean molecular weight of 
each component is imprinted in each dispersion relation. The difference between the Figure 
1 and 2 is that the growth rate at the larger scale decreases, compared with the original 
modes of Figure 1. This is because the suppression by the friction is more efficient at rather 
larger scale. The longer distance particles move, the more frequent the collision from which 
the friction originates is. It is noticed that the current model corresponds to the case that 
the fluid moves along the fleld line if there exists the magnetic fleld. Figure [2] shows that 
only the friction suppresses the instability along the fleld line, while the magnetic fleld does 
not. 

Figure [3] shows the growth rate when the friction is the same as Figure but there 
is the flnite magnetic fleld [i.e., vx/vsn = 0.6, i^nio/(^np^^sn) = 0.03 ]. For comparison, the 
two thin curves which correspond to the results of Figure [1] are also displayed, as well as 
Figure [2l A thick dashed curve is a mode, which is characterized by /xi. The reason why 
the curve is interpreted to relating to /Xi is that the critical wavelength is enlarged owing 
to the stabilization of the magnetic fleld [see the criterion (lB2p ]. That is, the magnetic 
fleld suppresses a kind of mode of the thermal instability. A thick dot-dashed curve is 
characterized by /^n- The critical wavelength of this mode is not affected by the magnetic 
fleld. Thus, the mode characterized by /in is hardly affected by the fleld. This occurs since 
the amount of the neutral component signiflcantly dominates that of the ion component. 

Figure H] shows the growth rate when the magnetic fleld is thousand times of that in 
Figure [3l and friction is ten times of that in Figure [3] [i.e., v^/vsn = 600, I'nio / {k^p^sn) = 
0.3]. The large magnetic fleld completely stabilizes a mode characterized by fii so that 
this mode does not appear in this flgure. Thus, there is displayed only one mode of the 
thermal instability characterized by /in. This is expressed by the thick dot-dashed curve. 
For comparison, the two thin curves are displayed as well as Figure [2l Figure [5] shows the 
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growth rate when the magnetic field is ten times of that in Figure H] and friction is the same 
as in Figure m [i.e., va/vsu = 6000, i^nio/ {f^npVsn) = 0.3]. The large magnetic field completely 
stabilizes a mode characterized by fii. The thick dot-dashed curve expresses the mode of 
thermal instability characterized by fi^i- For comparison, the two thin curves are displayed as 
well as Figure [2l It is noted that Figures H] and [5] show that the suppression by the magnetic 
field is saturated (see subsection 14.31) . 

In Figure El the friction is stronger than that in Figure HI but the magnetic field is the 
same [i.e., va/vsu = 600, ^'mo/l^np'^sn) = 5]. For comparison, the two thin curves are displayed 
like Figure [2l A mode characterized by /ii is completely stabilized, and then it disappears 
here. A mode characterized by /in is still unstable, while the growth rate is reduced. That is, 
even in weakly ionized fiuid, the growth of the thermal instability is affected by the magnetic 
field through the friction. This decrement of the growth rate is significant at longer scale. We 
also confirm that the characteristic wavelength at the maximum growth rate of the instability 
shifts shortward as the growth rate diminishes. This occurs since the friction is effective at 
larger scale and the conduction is effective at smaller scale, so that the growth rate at larger 
scale decreases significantly. Then, the characteristic wavelength at the maximum growth 
rate becomes small. 

It should be emphasized that all the critical wavelengths of the unstable mode charac- 
terized by /in are equivalent among Figures [1]- [61 The frictional terms in the right hand side 
of the dispersion relation (HM are independent of the critical wavelength for the condensation 
mode, as mentioned in subsection 13.11 The critical wavelength of the mode relating /in is 
not changed by the magnetic field with the friction as long as the thermal conduction is not 
changed by the magnetic field. This point is investigated more in the following subsection 

m 

4. DISCUSSION 

We study the basic property of the thermal instability of the weakly ionized plasma in 
the previous sections. In this section, basing on our understanding the thermal instability, 
we shall discuss the structure formation. 

4.1. Critical Wavelength and Thermal Conduction 

First in this section, the critical wavelength of the weakly ionized plasma in the magnetic 
field is discussed. As mentioned in section [3l the critical wavelength of the mode relating /in 
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is not changed by the magnetic field with the friction as long as the thermal conduction is 
not changed by the magnetic field. In fact, the thermal conductivity is varied owing to the 
magnetic field. In this meaning, we may oversimplify the problem in this paper. 

However, the thermal conductivity will decrease as the magnetic field strength increases, 
and then the critical length become smaller, according to equation (lB4p . Thermal conduction 
erases the temperature perturbation and thermally stabilizes the system, especially at smaller 
scale. The larger thermal cond uction make s the critical wavelength longer, owing to the 



criterion IB II or [B2| according to iFieldl (119651 ). On the other hand, the magnetic field lessens 
the thermal conduction of the ion component, because the ion component winds around a 
magnetic field. 

Therefore, we can insist that the critical length for the thermal instability of the weakly 
ionized plasma can not be enlarged by the friction and the magnetic field at least, while the 
effect of the magnetic field can infiuence not only the ion but the neutral component via the 
ion-neutral friction. This statement has universal validity for the thermal instability of the 
weakly ionized plasma. The effect of the friction and the magnetic field is mainly reduction 
of the growth rate of the thermal instability. On the other hand, the magnetic field directly 
enlarges the critical wavelength of the mode characterized by fii as expected in the one fiuid 
plasma model. 



4.2. Spatial Distribution of Magnetized Interstellar Medium 

In the previous section, we learn that the suppression of the thermal instability is char- 
acterized by the ion-neutral friction, strength of the magnetic field, and the direction of 
the field. Interestingly, the suppression of the instability is different owing to the magnetic 
field direction. This suggests the anisotropic nature of the suppression of the instability is 
imprinted in the morphology of the structure. In this standing point, we shall discuss how 
the structure formation proceeds in the magnetized ISM. 



4-2.1. Growth Rate Parallel to the Magnetic Field 

In an usual one fiuid approach, the motion of the plasma parallel to the magnetic field 
is regarded not to be affected. However, Figure [2] suggests that the growth rate of the 
thermal instability at larger scale decreases via the ion-neutral friction. This suppression 
of the instability yields the delay of the evolution in comparison to the evolution under the 
MHD approximation without the ion-neutral friction. 
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4-2.2. Distribution of Partially Ionized Plasma Perpendicular to the Magnetic Field Line 

Here, we suggest that the spatial distribution of the partially ionized plasma is elongated, 
whose semi-major axis is vertical to the magnetic field line. Figures [2] and [3l both of which 
friction strengths are the same, show the feasibility. The thermal instability in Figure [2] 
corresponds to the mode along the magnetic field line since B = 0, while that in Figure [3] is 
vertical to the magnetic field line. The growth of thermal instability vertical to the line is 
suppressed by both of the magnetic field and the friction, while the growth along the line is 
suppressed only by the friction. Thus, the growth rate along the magnetic field line is larger 
than that vertical to the line, and then the resultant morphology of the structure via the 
thermal instability is elongated. 

4.2.3. Difference of Spatial Distribution of Ion and Neutral Components 

We insist the growth rate depending on the direction of the field results in the morphol- 
ogy of interstellar structure. Furthermore, we learn in §3 the mode characterized by /i; are 
stabihzed by the friction and field, while that by /in can grow as found in Figure O Figure [H] 
shows that although the mode characterized by Hi is completely stabilized, the mode char- 
acterized by fin is still unstable. Consequently, we can insist that the neutral component 
condenses owing to the thermal instability, while the ion component still keeps spreading. 

It is noticed that the ion component may condense since it is dragged by the neutral 
component contracting owing to the thermal instability. However, once the distribution of 
the ion and the neutral components separates and the motion of the two becomes independent 
there, the ion component in the area is thermally stable owing to the magnetic field. This 
means the magnetic pressure lets the ion component be spread in the area even if the neutral 
component condenses. 

This thin ion component left behind by the contracting neutral component will be 
observed as a halo around H I clouds, and it emits, especially some forbidden lines. This halo 
possibly has wide distribution vertical to the magnetic field, as discussed in subsubsection 
14.2.21 This is because the ion component can condense, parallel to the magnetic field like the 
neutral. The magnetic pressure can not prevent the motion of the ion component efficiently 
as found in Figure [2l Thus, the morphology of the ion halo is also elongated according to 
the orientation of the field line. 
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4-2. 4- Origin of the Small dumpinesses of Partially Ionized Plasma 

It is suggested that there are small clumpy structures of ISM. Here, we examine whether 
the dumpiness originates via the thermal instability. According to Figures H] and O the 
growth rate of the thermal instability of partially ionized plasma depends on the strength 
of the friction. The scale length at the maximum growth rate of the instability becomes 
smaller, owing to the suppression via the friction which is effective at larger scale. Thus, 
we can expect the small dumpiness appears selectively via the thermal instability since the 
ion-neutral friction suppresses the large scale condensations. 

If we are interested in the small dumpiness of ISM, we should mention about the ther- 
mal conduction in relation to the magnetic field. According to subsection 14. H the small 
dumpiness may the results of the decrement of the conduction. It is noticed that the con- 
duction along the field line is not so affected, and then there is still ambiguity for the reduced 
conduction coefficient to perform the key role for the origin of the small dumpiness. 

It is also noticed that the thermal instability is a slow process relative to the dynamical 
processes. If molecular clouds form from the H I medium via the gravitational instability, 
the small dumpiness inside the molecular clouds appears after the formation of the clouds. 
This trend is enhanced if the effects of the magnetic field and friction are concerned. This is 
because the friction and the field suppress the large scale growth, while the short wavelength 
modes are still unstable. In a enough long span for the growth, the tiny structures can 
appear in ISM and affect the evolution of its structure especially as the non-linear effects. 
Thus, if we want to know how faint and small structures come to be observable, the growth 
time-scale of the thermal instability should be precisely examined. At this time, we never 
forget the effect of the magnetic field and the ion-neutral friction. 



4.3. Ineffectiveness of the Magnetic Field 

The modes in Figures H] and [5] have similar growth rate, while the magnetic fields are 
significantly different in each other. It implies that the suppression of the growth by increas- 
ing the magnetic field with the fixed friction has limits. By the way, a dispersion relation of 
only ionized fluid in a magnetic field, which is derived by equalizing the second bracket on 
the left-hand side in equation f|T^ to be zero, is reduced to 

+ u^VsikT = 0, (21) 

when k approaches 0. This is solved to be 



n = or n = — fgi^iT- 



(22) 
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Thus, we find that this solution is independent of a magnetic field. Even if the magnetic 
field becomes strong, the value of the growth rate at ~ is scarcely changed, where the 
friction is more effective than at small wavelength. Therefore, once the effect of the friction 
is fixed, increasing the magnetic field strength beyond a threshold level apparently becomes 
ineffective for the suppression of the growth. Not only information of the magnetic field but 
the friction is necessary to investigate the evolution of the thermal instability of the weakly 
ionized plasma. 

4.4. Comment on One Fluid Approximation 

The formulation of the fully one component plasma can not be always applicable to 
study the thermal instability of weekly ionized plasma. This is because there is the growing 
mode even when there exist the magnetic field and friction, whose growth rate is comparable 
to that of the neutral component when there are no magnetic field and friction. Indeed, this 
trend is found from the comparison among Figures [31 H] and O 

In addition, regarding weakly ionized plasma including the friction effect, the friction 
does not change the critical wavelength Ap. The friction tends to efficiently reduce the 
growth at much larger scale than the Ap. This feature yields that the critical wavelength in 
the weakly ionized plasma is not changed by the magnetic field via the friction, as mentioned 
in section [3] and subsection 14.11 Even if the magnetic field becomes strong, only the growth 
rate is reduced, keeping the critical wavelength constant. This is also the difference from 
the simple analysis of the one component magnetized plasma. The previous study of the 
one component plasma shows the stabilization at the smaller scales, while the current study 
insist Ap is unchanged even if the effect of the field is concerned. 

5. Application to Typical H I Region 

In this section, we shall study the origin of the small scale structure of H I cloud as 
an example. The thermal instability of the typical H I region with C ll cooling is briefly 
discussed here. 



5.1. Assumptions 

We assume that the temperature of H I is 100 K. The fluid mainly consists of H atom, 
and then fin = 1, /^j = 1/2- We assume that K = 4.9 x 10^ erg cm~^ s~^ K^^, consulting the 
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result of lUlmschneiderl (119701 ) as thermal conductivity. In addition, we adopt the collision 
cross- section as {av) = 2 x 10~^ cm^ s~^ for the typical friction value according to iNakano 
(Il976l ). We assume that interstellar pressure is about 10~^^ erg cm~^. Then, if the ISM 



is in equilibrium under the pressure of 10~^^ erg cm~^, the number density of the neutral 
component nneutrai and the magnetic field strength B are about 71.9 cm~^ and about 10~^ 
Gauss, respectively. The adiabatic 7 is assumed to be 5/3. 



According to IWolfire et al.l (119951 . Fig. 3 (b)), C ll cooling is dominant when T = 100 
K. Then we adopt the expression for the C ll cooling according to lWolfire et al.l (1l995l . their 
Appendix B). Assuming that the fluid consists of H atom, whose mass is expressed as mn, 
the cooling function A is 



A 



■ 2.54 X 10-^^Ac/cii ■ b^'^^HO + ■ exp 



92 
Y 



ergs s ^ g 



(23) 



Here, nx is the density of element X. 7^° and 7*^ are the coUisional de-excitation rate coef- 
ficients of C II with neutral hydrogen and electrons, respectively. The constant /en is the 
fraction of C I in C ll. Aq is the relative abundance of element C, d efined as nc, I n. n is the 
density of hydrogen nuclei, defined as n = nH+ + t^h"- According to lWolfire et al.l (119951 ) . 

HO 



7 



7 



2.1 X 10" 



T 
102 



-0.5 



8.86 X 10 



T 



-i°cm=^s^\ 



^0 + 0.484 — 
10^ 



+ 4.01 ( ^ 

10^ 



/ T 
3.39 — 

Vio^ 



(24) 
(25) 



and Ac = 3 X 10" 



The ionization degree and the fraction of C I in C ll, i.e., /en are ambiguous. Then, 
the ionization degree x and the constant /cn are hypothesized in some cases. We shall 
investigate the instability when /cn = 10"^ and /ch = 10~^. The more the amount of C 
II increases, the more efficient the cooling becomes. The io nization degree y i s assu med to 
be 10^^ according to the results around wh ~ 100 cm~'^ by iKamaya fc Nishil (120001 . Fig.l). 
We also investigate the instability on the extreme case when x = 10"^. In addition, we 
assume that A 



np 



Aip, AnT 



AiT, and K^q = Kiq. This is because the ionization degree 
is very small and the mixing of the neutral and the ion components is sufficient. Then the 
temperatures for the neutral and T\ for the ion become equivalent much faster than the 
cooling time. For example, if the temperature of the system changes by AT during the slight 
time At, the energy changes for the neutral and the ion are expressed as. 



and 



n^k-QAT ~ pnAnAt 



nifceAT ~ piAjAt. 
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From these equations, 

rin Ui 

Therefore, owing to mn ~ for H atoms, it is derived that An ~ Ai. 



5.2. Results and Implications 

Figure [7] shows the growth rate of the thermal instability when x = 10~^ and fcii = 10"^. 
The critical wavelength for the instability is estimated to be about 0.0135 pc owing to the 
criterion It is noted that = 0.920, = 920000, f3^ = 2.19 ■ 10"^ A = 2.73 • 10"^ 

= 219, = 0.775, = x/ fJ'ni = 5 ■ 10"'', according to our notation summarized in 
Appendix [Cl Figure [8] shows the growth rate when x = 10~^ and feu = 10~^. The critical 



wavelength for the instability is estimated to be about 0.135 pc owing to the criterion (IBip 



Here, = 0.920, = 920000, = 2.19 • 10-^ A = 2.73 • 10"^ = 219, = 77.5, 
Kp = x/ fJ'ni — ^ ' 10"''. Figure [3 shows the growth rate when x = 10~^ and fen = 10"^. 
The critical wavelength for the instability is estimated to be about 0.00239 pc owing to the 
criterion dE]). We adopt = 0.517, = 51.7, = 1.15■10-^ A = 1-44 -10"^ = 2.19, 
Cu = 1470, Kp = x//^ni = 5 ■ 10"'^. To understand more, in each set of these figures, the 
growth rates of the following three cases are also plotted: (1) the magnetic field is hundred 
times as the typical value, (2) the friction is hundred times as the typical value, and (3) 
both of the magnetic field and the friction are hundred times as the typical value. It is 
noticed that the value of the horizontal and the vertical axes is normalized by k^p, which is 
proportional to fcii- The more C ll (i.e., /cii) increases, the shorter the critical wavelength 
becomes and the more the growth rate increases. 

The each panel (b) in Figures [71 [H and [9] shows characteristic results. In these pan- 
els, we find some characteristic feature of the example. At larger wavenumber than fcmax, 
which corresponds to the maximum growth rate, the instability grows up as if there is no 
friction and field. We insist that the critical wavelengths of these cases are not enlarged 
because the neutral component does not suffer the effect of the magnetic field directly. At 
smaller wavenumber than /cmax, the growth rate is reduced significantly. That is, the growth 
of the instability in the significantly large-scale is suppressed. This may indicate the qua- 
sistatic evolution like the gravitational contraction via the so-called ambipolar diffusion, if 
the condensation due to the thermal instability is balanced by the magnetic field via the 
friction. 

Except the cases of x = 10^^, the overall trend of the results is the same as mentioned 
in section [31 The instability in the very low ionization degree cases of x = 10^^ is hardly 
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suppressed. In Figures [7] and [H the growth at only the large scale can be suppressed by the 
even large friction and magnetic field. In the case of x = 10"^ (Figure [9]), the ionization 
degree is higher than the other cases, then the ion- neutral friction is more efficient. The 
growth rate in the panel (d) of Figure M is reduced all over wavelength because of the large 
friction. The growth of the mode in the panel (e) in Figure [9] is extremely suppressed all over 
the wavelength because of the large friction and the magnetic field. Especially, the growth 
rate at the significantly large scale (i.e., near the origin of wavenumber) approaches to the 
zero owing to the friction. This figure shows the possibility that the thermal instability 
is suppressed as well as in Figure |6l The magnetic field can become strong for the ion 
component to be frozen in, and the friction between the ion and the neutral components can 
be strong even for the high density, which are expected at the final stage of the formation. 

It is noticed that under the conditions adopted in this section, the neutral component can 
always satisfy the thermal instability criterion, while the ion component does not satisfy the 
criterion. Importantly, the ionization degree is very small and then the ion component is too 
rare although the cooling owing to the ion component is effective for the thermal instability. 
Then, the thermal instability of the partially ionized plasma can not be suppressed except 
for the extreme case such as in the panel (e) of Figure O The partially ionized plasma can 
start to contract via the thermal instability in typical H I region, even when there exists a 
magnetic field. Even if the magnetic field becomes large, the growth rate is hardly suppressed 
except for the extreme case. 

6. SUMMARY 

The thermal instability of the weakly ionized fiuid is investigated with a linear pertur- 
bation analysis. The plasma is assumed to consist of two fluids of the ion and the neutral 
components. With this approach, the effect of the friction between the ion and the neu- 
tral components and the magnetic field are important. Here, the properties of the thermal 
instability of weakly ionized plasma and the observational implications are summarized. 

1. Modes relating to are not stabilized by the magnetic field. This is because the 
neutral component feels the field via the ion-neutral friction. 

2. Modes relating to /ii are directly stabilized by the magnetic field. This means that 
when the magnetic field is large, the ion component is hard to follow the motion of the 
neutrals which condensate owing to the thermal instability. 

3. The growth rate of the mode vertical to the magnetic field is reduced by the magnetic 
field and the ion-neutral friction. 
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4. The growth rate of the mode along the magnetic field decreases owing to the ion-neutral 
friction, especially at large scale. The suppression of the instability is more ineffective 
than that of the mode vertical to the magnetic field. 

5. The critical wavelength for the thermal instability is not affected by the friction. The 
effect owing to the two fiuid approximation is the reduction of the growth rate of the 
thermal instability. The magnetic field makes the critical wavelength of modes relating 
/li larger. 

6. To study the thermal instability of the partially ionized plasma, one fiuid approxima- 
tion is not always useful. 

7. The ion-neutral friction and the magnetic field affect the distribution or morphology 
of ISM, especially after the long time compared to the free-fall time. The difference 
between the growth rate along and perpendicular to the magnetic field is important. 
The partially ionized plasma possibly is elongated perpendicular to the magnetic field. 
In addition, the neutral component and the ion component of weakly ionized fiuid 
in the magnetic field are possibly separated each other. Then, the neutral compo- 
nent condenses owing to thermal instability, while the ion component left behind by 
the contracting neutral component still keeps spreading. Therefore, it is possibly ob- 
served that the weakly ionized fluid is elongated vertically to the magnetic field and 
surrounded by a halo which is rare and diffuse ions and emits forbidden lines. 

To conclude, the treating both of the independent motions of the neutral and the ion 
components yield the ion-neutral friction. The friction is important for the thermal in- 
stability of weakly ionized fiuid, especially when studying the structure formation such as 
molecular cloud at its final stage of formation. Our study indicates that the fully ionized 
plasma approximation or totally neutral fiuid approximation is not always applicable to 
weakly ionized plasma. 
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Appendix 

A. Derivation of the Dispersion Relation by Linear Perturbation Analysis 

We derive the dispersion relation (1141) in this Appendix. We adopt the perturbations as 
equation (|TT|) . that is, a(r, t) = ai exp{nt + ik-r). Basic equations ([I])-© are then linearized 
as follows: 

pnin + ipnoVni ■ k = 0, (Al) 

npnOVnl = -Pnlik - PnOl^niiVnl - Vil), (A2) 
Pnin '^—npnl = -pno(PnlAnp + TniAnT) -^n0^nlfc^ (A3) 



(A4) 



7-1 7 - 1 PnO 

Pnl _ Pnl T^l 

' Ttt 5 

PnO PnO nO 

pun + ipiov-a ■ k = 0, (A5) 

npiov-a = -piiik - Pnot'ni(f ii - *^ni) + -7-i{f^ X ^i) X Bo, (A6) 

Pun ^—npii = -pio(piiAip + TiiAiT) - KioTii/c^, (A7) 



7-1 7 - 1 Pio 

Pii Pii , Ti 



Pio Pio Ti 



- + ^, (A^ 



iO 



and 



nB, = t[vi,{Bo-k)-Bo{vi,-k)]. (A9) 



Using conditions Bq ± fc, || vh, k \\ v^i, and Bq \\ Bi, this multi-dimensional issue 
is reduced to one dimensional issue. Then, vanishing of the determinant of coefficients, a 
dispersion relation for partially ionized plasma is derived as follows: 
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PiO + PnO 47r PiO + PnO / 

k fsTi k. 
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A; A; 
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k 



nK 
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f Bq 1 , 2 , PiO + PnO , 2 



V 47r PiO + PnO PiO + PnO 



Vsnk^Vsik^ + Vsnk^VsikiT + Vsik-^V^nknT + t'snA:nT'ysiA;iT 

A;nK A;iK k^K kiK 



PiO 



PiO + PnO 



k'^Vsikip{VsnknT + 'f^snA^ 



-k'^VsnknpiVsikiT + Vgikj^ 



PiO + PnO 



k 

knK 
k 

kiK 



(AlO) 



We emphasize that this dispersion relation is not yet assumed to be weakly ionized 
plasma, such as pio <S Pno- Here are 



2 PnO 2 PiO 2 Bq 
Vsn = 7 , %i = 7 — , f^A = 1 , 

PnO PiO 47rpio 



Pn(7 - l)PnoA„p _ 
l^nT — 



k — 



Pn(7 - l)AnT , TlVsnPnO 
— , KnK — 



Pn(7 - l)^no' 



Pi(7 - l)PioAip , 



- 7^^,i^io 



Pi (7 - l)AiT , 'T^t'siPiO 
KiK — 



Pi(7 - i)-frio' 
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A -i^^A A -(^^'' 



Pn 



OA; 

dpi 



AiT = 



dAi 



pi 



and 



t'niO = luPiO- 



We assume that the fluid is weakly ionized (i.e., pio ^ Pno) and that temperatures of 
neutral and ion components are the same in the equilibrium (i.e., T^q = Tio). According 
to p = TZpT/p, assumptions pio <^ Pno and Tno = Tiq lead to pio <^ Pno- Using these 
approximations and the definition x = Pio/Pno, we can replace several expressions as follows: 



PiO + PnO 1 

— 

Pio X 

PnO 



PiO + PnO 7 
PiO 



4x 



PiO + PnO 7 

PiO + PnO 2 

7 ^^sn 

PiO + PnO 

PnO + Pi07 ^^sn 



PiO + PnO 7 

PiO + Pn07 2 
~ V 

PiO + PnO " 



and 



4vr PiO + PnO 



Then, the dispersion relation flAlOp is approximated to be 
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vlk^X + -Vsnk'^) 

1 A; 1 k 

X{'(JsikfVsikip{VsnknT + t^sn^:! ) {VsnkfVsnknpiVsikiT + V^ik- ) 

7 «^nK 7 «^iK . 



This is the dispersion relation 



B. Properties of the Dispersion Relation 



The pro perties of th ermal instability of one fluid approximation is briefly mentioned. 
According to iFieldl (119651 ) . the criterion of the condensation mode of the instability for the 
neutral component with or without the magnetic fleld is 



dK _ PnO dAn 



< 



nO 



PnO 



(Bl) 



If the cooling becomes eflicient with temperature diminishing, the system is thermally unsta- 
ble. This is expressed at the flrst term on the left-hand side of equation fIBll) . If the cooling 
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becomes efficient according to tlie contraction, tlie system is also tliermally unstable. Tliis 
is expressed at tlie second term on the left-hand side. The thermal conduction smoothes the 
perturbation to suppress the instability. It is efficient at small-scale events. This is expressed 
at the term on the right-hand side. 

Mathematically, whether the system is stable or not depends on whether the sign of the 
real part of the growth rate n is positive or negative. The system is unstable when the sign 
of the real part of n is positive. The sign of the real part of the solution of the dispersion 
relation is determined via algebra. 

Similarly, the criterion of the condensation mode of the instability of the ion in the 
magnetic field is 



Pio 9Ai 
^io dpi 



< -k 



PiO 



1 + 7- 



:b2) 



If dA[/dT[ > 0, the magnetic field always suppresses the instability. When dAi/dTi < 0, the 
criterion (lB2p is rewritten as 



^io dpi 



< 



PiO 



1+7 



(B3) 



Then, the effect of the magnetic field on the instability depends on the sign of dAi/dTi + 
k'^Kio/piQ. We note that this depends on the wave number k . 

The critical wavelength, i.e., Ap, above which the system is thermally unstable, is derived 
from equations (IBll) for the neutral, or (]B2p for the ion, respectively. The critical wavelength 
for the ion component depends on the magnetic field, while that for the neutral does not. 
The criterion (IBip is rewritten as 



Af 
27r 



1 
k 



KnO ( PnO <9An C^Ar 



PnO 



^nO C^Pn 



(B4) 



the criterion flB2l) is rewritten as 



Af 
27r 



1 
k 




1 + 7- 



respectively. 



PiO 5Ai _ 9Ai 
[Tio dp, dT, 



1 + 7- 



(B5) 



Moreover, the effect of the thermal conduction is mentioned. First, Observing to the 
conduction term V(-ft'nVTn) in equation ([3]) or the right-hand side of the criterion (IBll) . the 
conduction is more effective inversely proportionally to square of A. Then, the conduction 



- 26 - 



becomes less effective at large scale. Second, according to the criterion fIBip . the instability 
of the condensation mode prefers 

PnO <9An ^ dAn 



The preference allows us to approximate 




(B6) 



According to equation ([3]), KnoTno/ {XF{e)Y corresponds to the thermal conduction and PnoAn 
to the cooling rate. This equation (]B6I) indicates that the critical wavelength AF(e) is deter- 
mined at the point where the effect of the thermal conduction is equivalent to that of the 
cooling. The larger the thermal conduction is, the more the system stabilizes especially at 
smaller scale, and then the critical wavelength is enlarged. 



C. Normalization of the Dispersion Relation 

We define the following quantities: 



and 



kip 

k ' 

"Tip 



According to = 7(7^/p)T 



Vsn V Z^i^nO ' 



(CI) 



(C2) 



Assuming Tio ~ Ti 



Vsn V fJ'UnO V /^i 

Using these quantities, the dispersion relation (|T^ is normahzed as follows: 

+ + X^(3n) + YX^ + — (an + X^(3^ - 1] 

7 

X jr^ + y2 ^^^.^^c.. + fx^Ax^^ + F(/XniX)2 + Y{K,^Xf 

+ - ( yUni/tptti + /ini— X^ ) lilJ'mXy + 7(k„^X)^] - - (/ini^) Vni/t 



(C3) 
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(C4) 



We numerically solve this normalized dispersion relation flC4p in the main text. 
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Fig. 1. — Growth rates of each condensation mode of the ion component (the dot-dashed 
curve) and the neutral component (the dashed one), when B = and z/nio = 0. The dot- 
dashed curve corresponds to fully ionized fluid and dashed one to totally neutral fluid. The 
horizontal axis is the normalized wave number k/k^p. The vertical axis corresponds to the 
normalized growth rate ^[n/ {k^pVsn}]- 
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Fig. 2. — Growth rates of the condensation mode for the instabihty characterized by 
/in (the thick dashed curve) and /ii (the thick dot-dashed curve), when va/vsu = and 
^nio/(^np'ysn) = 0.03. The horizontal axis is the normahzed wave number k/k^p. The vertical 
axis corresponds to the normalized growth rate 3fJ[n/(fcnpfsn)]- This figure also corresponds 
to the case that the fluid moves along the magnetic field line. For comparison, there are 
also displayed condensation modes of the ion component (the thin dot-dashed curve) and 
the neutral component (the thin dashed curve) when B = and t'nio = 0, as well as Figure 
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Fig. 3. — Growth rates of the condensation mode for the instabihty characterized by 
/in (the thick dot-dashed curve) and /ii (the thick dashed curve), when va/vsu = 0.6 and 
^'nio/(^np'i^sn) = 0.03. The horizontal axis is the normahzed wave number k/k^p. The vertical 
axis corresponds to the normalized growth rate ^[n/ (fcnp^sn)]- The difference between Figures 
[2] and [3] is the magnetic field strength. For comparison, there are also displayed condensation 
modes of the ion component (the thin dot-dashed curve) and the neutral component (the 
thin dashed curve) when B = and z/nio = 0, like Figure 2. 
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Fig. 4. — Growth rates of the condensation mode for the instabihty characterized by /in 
(the thick dot-dashed curve), when vx/vga — 600 and / {kapVsr^ — 0.3. The mode of is 
completely stabilized by the magnetic field and does not appear in this figure. The horizontal 
axis is the normalized wave number k/k^p. The vertical axis corresponds to the normalized 
growth rate 3?[n/(A;npfsn)]- For comparison, there are also displayed condensation modes of 
the ion component (the thin dot-dashed curve) and the neutral component (the thin dashed 
curve) when 5 = and v-aio = like Figure 2. 
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Fig. 5. — Growth rates of the condensation mode for the instabihty characterized by /in 
(the thick dot-dashed curve), when va/vsu = 6000 and ^'nio/(^np'ysn) = 0.3. The mode of 
//i is completely stabilized by the magnetic field and does not appear in this figure. The 
horizontal axis is the normalized wave number k/k^p- The vertical axis corresponds to 
the normalized growth rate ^[n/{knpVsn)]- The difference between Figures H] and [5] is the 
magnetic field strength. For comparison, there are also displayed condensation modes of 
the ion component (the thin dot-dashed curve) and the neutral component (the thin dashed 
curve) when B = and z^nio = like Figure 2. 
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Fig. 6. — Growth rates of the condensation mode for the instabihty characterized by /in 
(the thick dot-dashed curve), when va/vsu = 600 and ^'nio/l^np'^sn) = 5. The mode of 
/ii is completely stabilized by the magnetic field and does not appear in this figure. The 
horizontal axis is the normalized wave number k/k^p. The vertical axis corresponds to the 
normalized growth rate 3fJ[n/(/Cnpfsn)]- The difference between Figure [Hand Figure [6] is the 
friction strength. For comparison, there are also displayed condensation modes of the ion 
component (the thin dot-dashed curve) and the neutral component (the thin dashed curve) 
when B = and u^iq = 0. 
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Fig. 7. — Growth rates of the condensation mode for the instabihty in the typical H I region, 
where T = 100 K and nneutrai = 71.9 cm~^. These panels show the case when feu = 10~^ and 
X = 10~^. The other parameters are assumed as (av) = 2 x 10~^ cm^ s~^, K = 4.9 x 10^ erg 
cm~^ s~^ K~^, Ac = 3 X 10~^, fin = 1, A^i = 1/2, 7 = 5/3. These panels are different in the 
strength of the magnetic field B and the friction, (a) both of B and the friction are zero; (b) 
both of B and the friction are the assuming typical values, i.e., B = 10~^ Gauss (the thick 
dot-dashed curve); (c) B is hundred times as the assuming typical value, while the friction 
the assuming typical value (the thick dot-dashed curve); (d) B is the typical value, while 
the friction is hundred times as the assuming typical value (the thick dot-dashed curve); (e) 
both of B and the friction are hundred times as the assuming typical values (the thick dot- 
dashed curve). For comparison, there are also plotted the condensation mode when B = 
and i^nio = (the thin dot-dashed curve) in each panel but (a). The horizontal axis is the 
normalized wave number k/knp- The vertical axis corresponds to the normalized growth rate 
^[n/{knpVsn)]- It is evaluated that k^p = 2.43 pc~^ and k^pV^n = 9.28 x 10~^^ s. It is noted 
that a„ = 0.920, = 920000, = 2.19 ■ 10-^ A = 2.73 ■ 10"^ k,,^ = 219, = 0.775, 
Kp = x/f^ni — 5 ■ 10"'', according to our notation summarized in Appendix O 
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Fig. 8. — Growth rates of the condensation mode for the instabihty in the typical H I region, 
where T = 100 K and rineutrai — 71.9 cm~^. These panels show the case when fen = 10^^ 
and X = 10^^. The other parameters are assumed as (av) = 2 x 10^^ cni"^ s^^, K = 4.9 x 10^ 
erg cm^-*^ s~^ K~^, Ac = 3 x 10^^, /x„ = 1, /Xj = 1/2, 7 = 5/3. These panels are different in 
the strength of the magnetic field B and the friction, (a) both of B and the friction are zero; 
(b) both B and the friction is the assuming typical values, i.e., B = 10~^ Gauss (the thick 
dot-dashed curve); (c) B is hundred times as the assuming typical value, while the friction 
the assuming typical value (the thick dot-dashed curve); (d) B is the typical value, while 
the friction is hundred times as the assuming typical value (the thick dot-dashed curve); 
(e) both of B and the friction are hundred times as the assuming typical values (the thick 
dot-dashed curve). For comparison, there are also plotted the condensation mode when 
B = and z/nio = (the thin dot-dashed curve) in each panel but (a). The horizontal 
axis is the normalized wave number k/kap. The vertical axis corresponds to the normalized 
growth rate ^[n/ {kr^pVsn)]- It is evaluated that k-^p = 0.0243 pc~^ and knpVgn = 9.28 x 10~^^ 
s. Here, = 0.920, = 920000, = 2.19 • 10-^ A = 2.73 • 10"^ = 219, = 77.5, 
= x//^ni = 5 • 10-^. 
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Fig. 9. — Growth rates of the condensation mode for the instabihty in the typical H I region, 
where T = 100 K and rineutrai — 71.9 cm^^. These panels show the case when fen = 10~^ and 
X = 10^^. The other parameters are assumed as (av) = 2 x 10^^ cm^ s^"*^, K = 4.9 x 10^ erg 
cm^^ s^^ K~^, Ac = 3 X 10^"*, fin = 1, A^i = 1/2) 7 = 5/3. These panels are different in the 
strength of the magnetic field B and the friction, (a) both of B and the friction are zero; (b) 
both of B and the friction are the assuming typical values, i.e., B — 10~^ Gauss (the thick 
dot-dashed curve); (c) B is hundred times as the assuming typical value, while the friction 
the assuming typical value (the thick dot-dashed curve); (d) B is the typical value, while 
the friction is hundred times as the assuming typical value (the thick dot-dashed curve); 
(e) both of B and the friction are hundred times as the assuming typical values (the thick 
dot-dashed curve). This case is the extreme one. For comparison, there are also plotted the 
condensation mode when B = and z/nio = (the thin dot-dashed curve) in each panel but 
(a). The horizontal axis is the normalized wave number k/k^^p. The vertical axis corresponds 
to the normalized growth rate ^[n/ {knpVsn}]- It is evaluated that knp = 12.8 pc~^ and 
knpVsn = 4.90 X 10-^=^ s. Wc adopt a„ = 0.517, a; = 51.7, = 1.15 • 10-^ A = 1-44 • 10"^ 
= 2.19, a = 1470, Kp = x//^ni = 5 • 10-^ 



